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We review recent studies on the measures of zero temperature quantum correlations 
namely, the quantum entanglement (concurrence) and discord present in the final state of 
a transverse XY spin chain following a quench through quantum critical points; the aim 
of these studies is to explore the scaling of the above quantities as a function of the quench 
rate. A comparative study between the concurrence and the quantum discord shows that 
their behavior is qualitatively the same though there are quantitative differences. For 
the present model, the scaling of both the quantities are given by the scaling of the 
density of the defect present in the final state though one can not find a closed form 
expression for the discord. We also extend our study of quantum discord to a transverse 
Ising chain in the presence of a three spin interaction. Finally, we present a study of 
the dynamical evolution of quantum discord and concurrence when two central qubits, 
initially prepared in a Werner state, are coupled to the environmental XY spin chain 
which is driven through quantum critical points. The qualitative behavior of quantum 
discord and concurrence are found to be similar as that of the decoherence factor. 

Keywords: quantum phase transitions; non-equilibrium dynamics; concurrence; discord. 



1. Introduction 

Quantum phase transition (QPT) is a zero-temperature transition of a quantum 
many body system driven by quantum flu ctuation s arising due to the presence of 
non-commuting terms in the Hamiltonian l l | 2 | 3 | 4 | QPT is generally associated 
with a diverging length scale (correlation length, £) and a diverging time scale 
(relaxation time, £ T ). Both the £ and £ r diverge in a power law fashion as the 
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quantum critical point (QCP), at A = 0, is approached by varying a parameter A 
of the Hamiltonian; £ ~ X~ v and £ T ~ £ z . The notion of universality demands that 
the correlation length exponent v and the dynamical exponent z depend on the 
dimensionality of the system, symmetry of the order parameter and the nature of 
the critical fixed point. In recent years, QPTs have been observed experimentally 
in a large number of systems, f or exa mple, in optical lattices a Mott insulator to 
superfluid transition is observed EM3. 

Recently, there has been a plethora of studies d irect ed to understanding the 
non-equilibrium dynamics of quantum critical systems These theoretic al stud ies 
are inspired by recent experiments on dynamic of ultracold atomic gases Q2LLU A 
natural question to ask that what happens when a quantum system is driven across 
a QCP by changing some parameter of the quantum Hamiltonian at a particular 
rate. It has been argued that even if the system is prepared in the ground state 
initially, there will be excitations (defects) in the final state reached following a 
quench across the QCP. This can be attributed to the existence of a diverging time 
scale in the vicinity of the QCP where the system becomes infinitely sluggish and 
hence can not evolve adiabatically however slow the variation of the parameter be! 
If a parameter A is varied as A = t/r, where r is the inverse rate of quenching, 
the defect density (n) in the final state is expected to satisfy a universal scaling 
relation given by h ~ T —vd/(vz+i) w j iere ^ j s the spatial dimension and v, z are 
the quantum critkal exponents as defined above; this is known as the Kibble-Zurek 
(KZ) scaling ^ ^ th at has been extensively studied for transverse XY spin chain 
in recent year j^- ^^ and has also been ex tended to different quenching scheme s, 
e.g., quenching through a multicritical point HSHSl (f or recent reviews, see ^^^ J). 

Ov er the past few years, the connection betw een between quantum informations 
and QPTs has been explored extensively It is interesting that quantum 
information theoretic measures can capture the ground state singularities associated 
with a QPT. The quantum correlations of a state can be quantified in terms of bi- 
partite entanglements. While the entanglement is a measure of the correlation based 
on the separability of two subsystems of a composite system, the quantum discord is 
based on the measurement on one of the subsyst ems. Both the concurrence (wh ich 
is one of the measures of bipartite entanglement) B 4 | 25 | 26 | an( ^ q uan t um discord 
show distinctive behaviors close to the QCP of a one-dimensional transverse XY 
spin model and also interesting scaling relations which incorporate the information 
about the universality of the associated QPT. For a spin chain i n the yicinit y of 
a QCP, behavior and scaling of both concurrence ^31 anc i discord I 28 | 29 | 30 131j f or 
example, a second order derivative of quantum discord with respect to the driven 
parameter of the quantum Hamiltonian shows a peak at the QCP. We note that the 
behavior of entanglement entropy and Re nyi entropy h ave been studied for trans- 
verse XY spin chain for both equilibrium I 32 | 33|34 J351 an( j j n ^e non-equilibrium 
state following a quench 

Here, we provide an introductory review of the recent studies which in fact 
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provide a bridge between the non-equilibrium dynamics of a quantum critical system 
and quantum information theory. One may, for example, raise the question what is 
the value of concurrence or discord in the final state of a quantum system following 
a quench across a QCP. If the dynamics is perfectly adiabatic, then no additional 
correlation is generated in the final state. However, as discussed above the passage 
through a QCP invariably leads to defects in the final state and these defects in turn 
lead to non-zero entanglement and discord ^Sl, i.e., non-adiabatic dynamics in 
the vicinity of the QCP results in non-zero quantum correlations. More importantly, 
both of them are found to scale with t in an identical fashion to that of the defect 
density. Finally, we show that equilibrium and non-equilibrium dynamics of a spin 
chai n can also alter the quantum correlation bet ween tw o qubits initially in a Werner 
state ^3 which are globally coupled to the chain I 40 | 41 |_ These observations open up 
a possibility of a deep underlying connection between non-equilibrium quantum 
critical dynamics and quantum information theory. 

The review is organized in the following way: in Sec. II we introduce the concept 
of quantum discord and also briefly that of concurrence. In Sec. III., we introduce 
the phase diagram of the transverse XY spin chain and transverse Ising spin chain 
with the three spin interaction and quenching schemes that are used in this review. 
In Sec. IV, the results for the transverse XY spin chain are presented and the 
possible scaling relation of concurrence and discord as a function of the quenching 
rate is analyzed; similar results for the three spin interacting model is given in the 
next section. In Sec. VI, we discuss the dynamical evolution of concurrence and 
discord between two qubits coupled to a transverse XY spin chain pointing out the 
similarity of their behavior with that of the decoherence factor. Concluding remarks 
are presented in Sec. 7. 

2. Quantum Discord and Concurrence 

The notion of quantum discord is based on the idea that a measurement perturbs 
the system and prepares the system in suitable eigenbasis of that particular mea- 
surement operator. Two classically equivalent expressions of mutual information 
become non-equivalent in the quantum case when one of them is associated with a 
local measurement. The difference between these two classically equivalent measures 
of mutual information is therefore the true indicator of quantum correlations present 
between two subsystems. This difference is in fact called the quantum discord. 

To illustrate this concept, let us start with a classical bipartite system composed 
of two subsystems A and B (see Fig. (fT])). Let us try to extract the information about 
the composite system AB from the overlapping part between H (A) and H(B). The 
mutual information associated with the composite system AB is quantified in terms 
of Shannon entropy H{p) where p is the probability distribution of the combined 
system. The classical mutual information is defined as 

I(p) = H(p A )+H(p B )-H(p), (1) 

where H(p 1 ),i = A,B stand for the entropy associated with the subsystem i with 
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H(A) H(B) 




H(A,B) 



Fig. 1. Venn diagram explaining the notion of mutual information. H(A) and H(B) represent 
information associated with subsystems A and B given in terms of the corresponding Shanon (von 
Neumann) entropies, respectively. H(A, B) represents information associated with the composite 
system AB. H(A\B) and H(B\A) denote information associated with A for a given information 
of B and vice-versa, respectively. H(A : B) signifies the mutual information associated with both 
A and B. 



probability p l . Alternatively, this mutual information can be denned through the 
relation 

J(p) = H(p A )-H(p\p B ), (2) 

where H(p\p B ) — H(p) — H{p B ) is the conditional entropy. In the classical situation, 
it is obvious that I — J as, any measurement on B leaves the state of the composite 
system unchanged. 

Let us now attempt to generalize this to the quantum case when the classical 
Shannon entropy gets replaced by the quantum von Neumann entropy expressed in 
terms of the density matrix. The natural quantum extension of Eq. ([1]) is 

I(p) = s(p A ) + s(p B ) s(p), (3) 

where p A and p B are the density matrices for the subsystems A and B, respectively, 
while p is that of the composite system. In the quantum context, the conditional en- 
tropy is defined through a local measurement carried over the subsystem B. Usually 
one employs a von Neumann type of measurement consisting of a set of one dimen- 
sional projection operators Following a local measurement on the subsystem 
B, the probability of the final state pk of the composite system where 

Pk = — (I ® B k )p(I ® Bk), (4) 
Pk 

is given by pk = tr(7 <£> Bk)p(I <8> -Bfc), where I is the identity operator defined 
in the subspace of A. One can therefore define the quantum generalization of the 
conditional entropy given by s(p\{Bk}) — J2kP kS (P k ) leading to the equivalent 
definition of the mutual information given by J(p\{Bk}) = s(p A ) — s(p\{Bk})- 
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The classical correlation can be ext rac ted by maximizing the above expression 
of expression of J in the following way S2l 

C(p)=max {I§k} J(p\{B k }). (5) 

As introduced by Olliver and Zurek 127-1, the difference between quantum mutual in- 
formation I(p) (Eq.([3])) and the measurement induced classical correlation (Eq.([5])). 
i.e., 

Q{p) = I(p) - C( P ) (6) 

is the quantum discord which is the true measures the quantum correlations. 

It is worth noting that I represents the total information (correlation) whereas 
C is the information of A following a measurement in B. If Q = 0, one can con- 
clude that the measurement has extracted all the information about the correlation 
between A and B i.e., system has no quantum correlation. On the other hand, a 
non-zero Q implies that entire information about the A can not be extracted by 
local measurement on B. This difference originates from the fact that the subsystem 
B gets disturbed in the process of measurement which does not happen in classical 
information theory. 

Let us quickly recap the notion of concurrence which is a measure of bipartite 
entanglement based on the separability approach. If states of the composite system 
could be written in terms of tensor product of states of the subsystems A and B then 
the concurrence is zero i.e., they are unentangled. The concurrence for a mixed state 
of two qubits (^4 and B) is given by C nc — maxjO, \fX\ — x/A^ — \Aa — a/XPK where 
Ai's are the eigenvalues of p(a y ® a y p*a y ®a y ) in decreasing order E4| 25 | 26| There 
could exist quantum states with zero entanglement but nonzero quantum discord. 
As to be shown below, the final state of the spin chain following a quantum quench 
may represent this type of a state if the rate of quenching is below a threshold value. 



3. Model and Quenching Scheme 
3.1. Transverse XY spin chain 

The Hamiltonian of one-dimensional spin-1/2 XY model in a transverse field with 
nearest neighbor ferromagnetic interactions is given by S3 

H = -\ £ [(1 + 7)4<4 +1 + (1 - lK° y +1 + hall (7) 

i 

where cr's are the Pauli spin matrices and the subscripts stand for the spin direction 
and superscripts are the lattice indices. The parameter h is the transverse magnetic 
field and 7 measures the anisotropic nature of the interactions; 7=1 corresponds 
to the transverse Ising chairM 

The model ([7]) can be exactly solved by mapping t he s pins to spinless non- 
interacting fermions via a Jordan- Wigner transformation SHI followed by a Bogoli- 
ubov transformation. The effective two-level Hamiltonian for each pair of momenta 
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Fig. 2. (Color online) The phase diagram of one dimensional XY model in a transverse field given 
by Hamiltonian The vertical red lines at h = ±1 denotes Ising transition from ferromagnetic 
phase to paramagnetic (PM) phase. The horizontal blue line stands for anisotropic phase transition 
between two ferromagnetic phases FM^ and FM S with magnetic ordering in x and y directions, 
respectively. 



±k in terms of the states |0) and \k, —k) = ctc_ JO) is given by 

_fh + cosk 7sinfc \ 
k ~\ j sink -{h + cosk) ) ' [ ' 

where c/.(c|,) is the Fourier transform of Ci(c\) — n}=-oo °z( — l) l(7 - |0) , and 
—k) denote the state with no and a pair of c-fermions, respectively. In the reduced 
Hilbert space, any general state can be represented as a superposition of |0) and 
\k, —k) with time dependent amplitudes Uk{t) and Vk(t) such that ipk(t) — Uk(t)\0) + 
Vk(t)\k, — k) . The phase diagram for the model is shown in Fig. The dynamical 
critical exponents and the correlation length exponents associated with the Ising 
transitions (at h — ±1) and the anisotropic transition (7 = 0, —1 < h < +1) are 
given by v = z = 1. 

We shall study the behavior of quantum discord in the final state after quenching 
the system across Ising crit ical po ints following the quenching scheme h(t) = t/r, 
with t going from — 00 to oo ^^ l The diverging relaxation time close to the QCPs 
at h = ±1 leads to defects in the final state. Initially (at t —> —00), the system is at 
the ground state |0) where all the spins are aligned in the —z direction. At t — > 00, 
on the other hand, the system is in an excited state due to non-adiabatic excitations 
in the vicinity of the QCPs at h = ±1; the probabilities of non adiabatic transition 
for the mode k is given by the Landau-Zener transition formula 

Pk = \uk(+oo)\ 2 = exp(-7TT7 2 sin 2 k). (9) 



In the limit r — > 00, only the modes close to the critical modes (k = or k = n) 
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Fig. 3. Equilibrium phase diagram of the three spin interacting Ising model {12) . Solid lines show 
phase boundaries between PM and FM phases. 



contribute to (|9|) and one arrives at a simplified form p^ — exp(— Trj 2 k 2 T). 

Wc further use a quenching scheme that drives the system across a qua ntum 
multicritical point (MCP) "A" . For a linear path, h and 7 are related as^^ 

fc(7) = l + |7(t)|asn(i); *y(t) = ~, (10) 

T 

with t varying from — 00 to +00, where the MCP is at t = 0. We investigate the 
scaling of discord in the final state following this quench using th e expr ession for 
the probability of excitation pk which in this case gets modified to 

Pk = exp(— 7tt(1 + cos/c) 2 sin 2 k). (11) 



3.2. Transverse Ising model with a three spin interaction 

Let us now introduce a one-dimen sional three spin interacting transverse Ising sys- 
tem described by the Hamiltonian 

i i 

where J x is the strength of the nearest neighbor ferromagnetic interaction and J3 
denotes the strength of the three spin interaction. In the limit J3 — > 0, the model 
reduces to transverse Ising model. Moreover, Hamiltonian (|12p can be mapped to a 
transverse XY spin chain with competing (ferro antiferromagnetic) inte ractions in 
the x and y components of the spin using a duality transformation SSI Therefore, 
it can be exactly solved using JW transformations in spite of the presence of the 
three spin interactions. 

The reduced Hamiltonian in the Fourier space is of the form 

(h(t) + J x cosk — J3COS2&;) sin fc — J3sin2A:) 

—i(J x sin k — J3 sin2fc) —(h(t) + J x cos k — J3 cos 2k) 
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We can analyze the gap of the spectrum and it can be seen that gap vanishes at 
h = J3 + I, and also at h = J3 — 1, with ordering wave vectors tt and 0, respectively; 
these transitions are equivalent to the Ising transitions in the transverse XY spin 
chain while the phase transition at h — — J 3 (and J3 > 0.5) belongs to the 
universality class of the anisotropic transition. The equilibrium phase diagram is 
shown in the Fig. © S^. 

Let us consider the quenching dynamics of the system across Ising critical points 
following the quench scheme h(t) — t/r with t going from — 00 to 00 The 
probability of non-adiabatic transition close to the qua ntum critical points are again 
given by the corresponding Landau-Zener formula 



4. Pairwise Correlations, Quantum Discord, and Concurrence 

Let us now address the question what would be concurrence and discord in the final 
paramagnetic phase of the Hamiltonians ([7]) following quantum quenches using the 
scheme h{t) = t/r. The initial state (t — > —00) is a direct product state with 
all the spins pointing in — z direction. If the dynamics is perfectly adiabatic, the 
final state is also expected to be a direct product state (with all spins pointing 
up) . This implies that there is no quantum correlation generated if the dynamics is 
perfectly adiabatic which is the never the case if the system crosses a QCP during 
the temporal evolution. Passage through a QCP generates defects in the final state 
which in turn leads to non-zero quantum correlations that show up in the measures 
like concurrence and discord. 

To calculate these measures, we use the two-spin density matrix for spins at the 
sites i and j — i + n. Let us recall the generic form of the density matrix given by 



where c\ = (<7*erj), C2 = (o^o^), C3 = {a\al), C4 = c§ = (er!,). For the Hamiltonian 
([7]), the equivalence of X and Y directions (in the spin space) demands c\ = c 2 . 
The density matrix can also be rewritten in the form 



Pk = exp[— 7rr(sin k — J3sin2/c) 2 ]. 



(13) 



1301471 




/ a% &y\ 
a% b% 
b%* a% 

\&r a n _) 



(15) 



where the matrix elements are given in terms of the two-spin correlation functions 
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in the following manner: 



al = -<(1 ± al)(l ± a\ +n )) = 1 + c 3 ± 2c 4 , 
«o = J<(l±^)(lT^+")) = l-c 3 , 



(16) 



Moreover, the up-down symmetry of the Hamiltonian leads to a simplified form of 
the density matrix whe re som e elements vanish S3. 
Defining a quantity 



0n 



one gets 



/ — p k cos(nk), 
Jo T 

c 4 = c 5 = (4) = 1 - 2/3 , 



4^ 



(17) 



(18) 



The expressions for ci and C2 can be computed for different value of n which we 
present below for n < 6: 

'^(l-2/3 ), 



ci = c 2 



(1 - 2/3 ) 2 /3 2 2 - 4/3f + f (1 - 2/3 ) 3 - 2/3|/3 4 (l - 2/? ), 



±[/3 6 {(l - 2/3 ) 2 - 4/3 2 )} + 4 /3 2 {/3 2 + /3| - &(1 - 2/3 )}] x [16/3 2 /3 4 + 
[ (1 - 2A)){(1 - 2/3 ) 2 - 8/3 2 - 4/3 2 }], n = 6 

The eigen values of the density matrix are obtained in terms of the correlators CiS 
as™ 



A = 


> 


+ c 3 ) + 


\/ 4c i + ( c i 


-c 2 ) 2 ] 


Ai = 




+ c 3 )- 


V 4c l + ( c l 


-c 2 ) 2 ] 


A 2 = 




-c 3 ) + 


(ci + c 2 )], 


and 


As = 




-c 3 )- 


(ci +c 2 )], 





(19) 



which can be expressed in terms of /3's using the equations ([To]) . (fT7|> and JT? 
Using Eq. (fT7|) we note that /3„ = for odd n as is invariant under k —> n — k. 
At the same time, (a±a±~ n ) = 6™ = for all n since the expectation values of 
a pair of fermionic annihilation or creation operators do always vanish. Moreover, 
( a ± a ^ n ) = ^ 2 = for odd n since the quantities 6 2 change sign under the Z 2 
transformation | lt) | 37 | 47 | 

On the other hand, 6 2 = C\ +c 2 for even n. The variation of mutual Information 
/, the classical correlation C and the quantum discord Q = I — C with t are 
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therefore studied for both critical and multicritical quenches (|10|) for even n. Let us 
assume that the spin at site i to be subsystem A and at site j as subsystem B. The 
reduced density matrix for the subsystems A and B can be expressed as 

1 



PA 

PB 



(P ® V + c 4 -T ®(J 3 Z ), and 

{P®P +aai ®P), 



(20) 



with eigenvalues 



A 4 = -(1 + c 4 ), and 
A 5 = ^(1 - c 4 ). 



(21) 



The total mutual information I(p) is expressed terms of von Neumann entropies, 
which when substituted in Eq. ([1]) gives 



Hp) = s (pa) + s(p B ) - ^2 Xa log2 Ac 



(22) 



where s(pa) — s(pb) = — A 4 log 2 A 4 — A5 log 2 A5. In order to calculate the classical 
correlation, we employ a set of projector for local measurement on the subsystem B 
given by Bk = VHkV 1 where LTfc = \k){k\ : k = +, — is the set of projectors on the 
computational basis |+) = ^(|0) + |1)), |-) = ^(|0) - |1)) and V € U(2) where 
V is parametrized over a Bloch sphere given by 



sinfe^ -cosf 



(23) 



where the polar angle 9 lies between and tt and the azimuthal angle <fi can be varied 



from to 2-7T. Following the standard techniques 
correlation by maximizing 



30 



C{p) = s{p A ) - s(p+), 



we can obtain the classical 
(24) 

-f). Below we summarize the 



where p+ is the density matrix for the outcome | fc) = 
final results for n = 2. 

The exact expressions for mutual information and classical correlation are given 
below. 

/ = -2(1 - A>)iog a (i - ft) + ((1 - M 2 - /3 2 2 )iog 2 ((i - M 2 - Pt) - 2A)fog 2 (A)) 

+ (/3 2 - (3 2 2 )logM - pi) + i{4/3 (l - h) + 4/31 + A(l - 2/?o)} 



x log 2 
x log 2 



{4/3 (l - /3 ) + 4/3| + /? 2 (1 - 2/3 )} 
{4/3 (l - /3 ) + 4/3| - /? 2 (1 - 2/3 )} 



+ -{Ap (l - A)) + 4/3 2 2 - - 2/3 )} 



(25) 
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and 



C = -(1 - A>) log 2 (l - fo) - p logaCSo) + \ ( 1 - (1 - 2/3 )y 1 + ^ 



x log 



x log 



II. /. /3 2 2 



, ,i-(i-2A,)yi+ , 



i+(i-2A))yi . 4 



i /. _„ . /. pi 



+ 2 U + (l-2A)Vl+ , 



(26) 



Similarly one can obtain the expressions for n = 4 and n = 6 using the appropriate 
equations. We note that I and C and hence Q = I — C depends entirely on /3's 
which are in turn dependent on the quench rate r _1 through the defect density pk- 
The spin chain has non zero concurrence! '^^ if |&Jg | > ^/°+ a -- The concurrence 
is given by C nc = max {0,2(\b%\- ^/a^)} . 



4.1. Comparative study between discord and concurrence 

We would like to make a comparative study between the discord and concurrence 
as a function of quenching rate r. The variation of quantum discord Q with r for 
n = 2, 4, 6 following a quench across the Ising critical point is shown in Fig. (|4"|. Q 
vanishes in both the limits r — ► (sudden limit) and r — > oo (adiabatic limit) due 
to the fact that the final state is nearly a direct product state in either cases. In the 
adiabatic limit, the system nearly reaches the expected final state with minimum 
defect while in the sudden limit, the final state is almost identical to the initial 
ground state (which is a direct product state). Therefore Q increases monotonically 
with r and reaches a peak value at an intermediate r = r m , and starts decreasing 
for t > T m . A similar behavior is also seen for von-Neumann entropy density as 
shown inl l6 | 17 | 

As the lattice spacing n increases, r m is also shifted to higher values and the mag- 
nitude of Q decreases indicating that the correlation generated through quenching 
is short-ranged. Classical correlation as defined in Eq. (|26|) also exhibits a qualita- 
tively identical behavior though it is smaller in magnitude in comparison to discord 
(Fig. {!])). It is worth noting, the classical correlations however show some ini- 
tial fluctuating behavior for r — > followed by the monotonic increase (see inset 
Fig. flU). This is due to the fact that the classical correlation becomes zero at 
some intermediate times when the measurement basis matches with the state of a 
subsystem, namely i-th spin. The value of C is found to be one order of magnitude 
less than Q implying that correlations present in the system are essentially quantum 
mechanical. 

The variation of concurrence in the final s tate for a quench across an Ising crit- 
ical point with r has been studied recently and the comparison is shown in 
Fig. (|4|. Although the variation of discord and concurrence are qualitatively simi- 
lar, we emphasize following differences. The magnitude of discord is less than that 
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Fig. 4. The variation of quantum discord Q with r for n = 2 (solid line), 4 (dotted line) and 

6 (dashed line) in the final state following a linear quench across the Ising critical points where 

7 = 1. Inset shows the variation of concurrence (C nc ) for same parameter value. 

0.0025 
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11.11005 



50 100 150 200 250 

X 

Fig. 5. (color online) Variation of classical correlation with r for n = 2 (solid line), 4 (dotted line) 
and 6(dashed line). Inset shows small peaks for r — > 0. 

of concurrence for the same value of n nearly by one order suggesting that discord 
represents quantum correlations. Moreover, Q shows a peak at a value of r m which 
is very small in comparison to the corresponding r m for concurrence. Based on the 
above observations, we conclude that the measurement based approach employed 
in calculating discord provides a quantitatively different result for quantum correla- 
tions. Moreover, the study on concurrence indicates the existence of a threshold 
value of t above which the bipartite entanglement is generated through quenching. 
On the contrary, discord is non-zero for all r as we observe negligible shift close to 
r = for different n as shown in Fig. ((4]). Therefore, in some situations quenching 
leads to a final state that happens an ideal example of a quantum state with zero 
concurrence but non-zero discord. 

4.2. Scaling of discord and concurrence 

We would now like to explore whether Q satisfies a universal scaling relation as a 
function of r; from the KZ scaling relation we recall that the defect density scales 
as \j\fr for quenching through the Ising critical points. It has been reported that 
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logQ 




2.8 2.85 2.9 2.95 3 3.05 3.1 3.15 3.2 3.25 3.3 3.35 

logx 



Fig. 6. The Variation of quantum discords Q as a function of t for n = 2 (solid line), 4 (dotted 
line) and 6 (dashed line) are shown on a log-scale following a quench across the Ising critical points 
with 7 = 1; the slope is 0.5 in all cases indicating that discord satisfies identical scaling relation 
as that defect density. Inset shows the variation of concurrence (C nc ) for n = 2 with the same 
parameter values ; the slope is again ^ 0.5 




-2 - 



-2.1 1 1 1 1 1 1 1 1 

2.65 2.7 2.75 2.8 2.85 2.9 2.95 3 

logx 

Fig. 7. The variation of discord for n = 2 following a multicritical quench along a linear path 
with slope —0.19 which matches well with the value of the exponent —1/6 obtained for the 
defect density . In the inset, we show the similar variation of concurrence with r for n = 2 and 
the slope is —0.13 which is in close agreement with the exponent —1/6. 

concurrence also scales as l/y/r in the limit of large -r ET] j n p^g (Jg]) we analyze the 
scaling of discord which clearly shows a 1/ ^fr fall. In order to explore the scaling of 
Q analytically in the limit r — > co, we analyze the asymptotic behavior of the terms 
of / (Eq. (|25jl ). The first two terms together decay as 1/r, while the contribution 
from the remaining four terms has been found to scale numerically as 1/y/r, which 
determine the scaling when r — > oo. Although, a closed power-law form is not 
obtained, our studies apparently points to the fact that discord does also satisfy a 
scaling analogous to that of concurrence or defect density. 

We further verify this claim by investigating the scaling of discord following 
a linear quench ing acr oss the MCP 'A' (see Eq. (9)) for which the defect density 
scales as r -1 / 6 . I 18 | 19 | j n pjg^ we present the scaling of discord (and also that 
of concurrence) with respect to r; our numerical results again indicates that the 
scaling of discord is likely to be same as that of the defect density. 
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Fig. 8. (color online)The Variation of quantum discords Q for n = 2 as a function of r with 
different values of J3. Discord decays slowly around J3 = 0.5 than any other values of J3. Inset 
shows the Variation of concurrence as a function of r with different values of J3 for the same 
lattice spacing n = 2. Concurrence becomes zero when J3 > 0.5. 
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Fig. 9. (color online)The Variation of Q for n = 2 as a function of J3 with different values of t. 
A sharp peak of discord gradually approaches around J3 = 0.5 as r increases. The inset shows the 
Variation of concurrence as a function of J3 for the above parameter values. Concurrence shows a 
peak around J3 = 0.5 and becomes zero as J3 > 0.5 for any value of t. 



5. Some Results for Transverse Ising spin chain with three Spin 
Interaction 

In this section, we extend our previous study (in Sec. IV.) to the context of a three 
spin interacting transverse Ising chain (|12p for quenching scheme h(t) = t/r. Using 
reduced matrix (| 1 5jl for the present model, one obtains the variation of quantum 
discord Q as a function of t with different values of J3 as shown in Fig. ©. The peak 
hight of discord decreases with increasing three spin interaction strength. This can 
be explained by noting that the magnitude of the defect density is inversely related 
to JfJ-^ and discord, as discussed previously, is expected to be proportional to the 
defect density. The Inset of Fig. ([8]) shows the variation of concurrence as a function 
of quench rate. Although the behavior of both of them are similar, the concurrence 
vanishes if J3 exceeds 0.5. One can further check that the discord and concurrence 
satisfy the same scaling as the defect density for both critical and multicritical 
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quenches. 

Fig. ([9]) shows the variation of quantum discord and concurrence as a function 
of J3 for a fixed value of quench time r. Discord shows a peak at a value of J3 that 
approaches J3 = 0.5 and becomes sharper with increasing r; on the other hand, 
concurrences vanishes for J3 > 0.5 for all values of r. 

6. Dynamical Evolution of Quantum Discord and Concurrence 

In the earlier section, we presented a review of the scaling of concurrence and discord 
in the final state following a quench across a QCP and their connection to the scaling 
of the defect density. In this section, on the other hand, we study the dynamical 
evolution of these measures between a pair of qubits which are coupled to a driven 
quantum spin chain; we choose the environmental spin chain to be the transverse 
XY spin chain defined in Eq. . This environmental spin chain is globally coupled 
to the two qubits (denoted by A and B) through the interaction Hamiltonian given 
by Hse — Yli + Sf)a \ where S denotes the coupling strength between 

system and environment | 40 | 41 ] ^ e note tnat ^gA _|_ S B ,a^] = (a — x,y,z) so 
that 5/2(S A + Sf) is conserved during time evolution. The global Hamiltonian of 
the system and qubits can then be written as, 

4 

H = Y,M(<P»\®H%r(h(t)+e^, (27) 

with I^H/z = 1, 2, 3, 4; 1^) = I + +), |p 2 ) = |--), |^,4> - 1/V2(\ + -)±\- +») 
denoting the pth eigenstate of the operator (S A + Sf), i.e., central two qubits 
system, with the p,th. eigenvalue e M (ei = +5, £2 = — <5, £3,4 = 0). The parameters 
hfi are given by h^(t) = h(t) + £ M , and the environmental Hamiltonian H E ^ is 
obtained from He by replacing h(t) with hpit). We therefore find that the coupling 
between the qubits and the environment, splits the Hamiltonian into two branches 
with transverse field (h(t) + S) and (h(t) — S), respectively. For our studies, employ 
the following time variation of of the transverse field h(t) given by hit) = 1 — t/r 
while t varies from —00 to +00. 

Let us also assume that the central qubit system and the environment (spin 
chain) are initially uncorrelated with the total density operator of the compos- 
ite system at t = —00, is given by p tot (—oo) — p AB (~oo) ® p E (—oo), where 
p AB {— 00) is the initial density operator of the two central qubits and \^e{— 00)) 
is the initial state of the environment yielding the density operator p E (— 00) = 
\tJje(— 00)) (iPe{— 00) \. The time evolution of the composite system is determined by 
the time evolution operator U(t), as p tot {t) = U(t)p tot (—oo)U(t)^ . The operator 
U (t) l^QI, generated by the Hamiltonian (|27|) is 

4 

t^HEi^x^i^m (28) 

u=l 
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with £7 B M (t) = exp(— iHJ 1 (t)). Finally, the reduced density matrix of the two central 
qubits at an instant t is obtained by tracing over the environmental degrees of 
freedom, i.e., p AB it) = Tr E [p tot {t)}. We further assume that the initial state of the 
two central qubits is a Werner state given by 

/*(-cx>) = ^I^ + a|<^|, (29) 



where \<f>) = (|++) + | ))/v / 2- Clearly, the werner state (|29p becomes totally 

mixed for a — 0, and reduces to a pure state \<j>) in the case of a = 1. 

Noting that the coupling between the qubits and the environment provides two 
channels of evolution of the environmental wave f unc tion, the reduced density matrix 
of the two central qubits can now be written as SH 



AB 



1 + a 





2ay/\D(T)\\ 
1 - a 
1 - a 



(30) 



\2ay/\D(t)\ l + a J 



where where the decoherence factor D(t) = \(ip+(t)\i/>-(t))\ 2 with \ifi + (t)) and 
|<0_(t)) being the states of environmental spin chain at tim e t evolving with the 
Hamiltonian Hgihit) + 5) and HE(h(t) — S), respectively SSI. 

To evaluate D(t), we rewrite the Hamiltonian (0 with modified h (due to the 
coupling S) in terms of JW fermions which then can be decoupled into a sum of 
independent (2 x 2) Hamiltonians in the Fourier space S31 j n t ne basis |0) and 
\k, —k), which represent no quasiparticle, and quasiparticles with momentum k and 
—k, respectively, the environmental Hamiltonian can be split as Hjj(t) = Ylk ^fc 
where 

7 sin k 
-(h(t) ±5 + cos k) _ 

The general wave function i/j±(t) can therefore be written as (see Sec l3.ll) 

hk(*)> =ni^* (*)> = n [«*(*)io)+^(t)iA,-*)] . (si) 

k k>0 

The coefficients and are obtained by solving the Schrodinger equation 
id I dt (u^ (i) , (t)) T = Hj*r(t) (uf(t),v^(t)) T where A T represents the trans- 
pose operation of the row matrix A. Hence, the expression of D{t) is given by 
U k Fk(t) = Uk KMHt) + S)\MHt) - <5))| 2 , or, 



h(t) ± S + cos k 
7 sin k 



D(t) = cxp 



N 
2^ 



dk InF, 



(32) 



where can be written in terms of u k and . One finds the exact form of Fk to 
be given by HHHS1 

-2nrk 2 -inrk 2 



F k (t) w l-4shr(4M) 



(33) 
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In weak coupling limit {6 —¥ 0), one can also get a closed form ex pressio n of deco- 
herence factor ([32]) after crossing the first QCP at h = 1 given by | 48 | 49 | 

( 8(V2-1)A^ 2 V 
D(t) « cxp ^- {V n J_ J C, (34) 

where £ is the contribution due to the fidelity factor for adiabatically evolving 
modes. 

0.6 



0.5 



0.4 



0.3 



0.2 
0.1 


-2 -1.5 -1 -0.5 0,5 1 1.5 2 2.5 3 

hft) 

Fig. 10. Quantum discord is plotted as a function of transverse field h{t) with three different 
values of a; a = 0.9(dashed line) corresponds to the upper one, a = 0.5(dotted line) corresponds to 
the middle one, and a = 0.3(solid line) is represented by the lower one. We have chosen N = 500, 
8 = 0.01 and r = 250. Inset shows the variation of concurrence as a function of transverse field 
h(t). 




We shall study below the behavior of discord and concurrence using different 
quenching schemes a function of a and h(t), as the environment is driven across 
QCPs. This study reveals that the quantum discord and concurrence follow qual- 
itatively similar kind of behavior as that of the decoherence factor of the central 
qubits. We show that the purity of the central spin state is not entirely determined 
by a, rather it does also depends on the decoherence factor. The central qubits will 
be in a pure state only when a — 1 as well as D(t) = 1. 

U sing the reduced matrix Eq. (|30[) . one can find the quantum discord in the 
form 

4 

Q[p AB {t)\ = — ^ log a (-^) + E A ™ l0 S2 Am + f(a), (35) 

m 

where A m 's are the eigen value of p AB {t) and /(a) is a function of a. Similarly, the 
concurrcnc can be obtained as 

C nc {p AB (t)}=max[a(^]DW\ + l)-l,0}. (36) 

The equilibrium behavior of discord between two central qubits coupled to a 
transverse spin chain has been extensively studied SB. We would like to investigate 
large S behavior of discord and concurrence numerically studying Eqs. (|35[) and (|36[) 
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Fig. 11. Quantum discord is plotted as a function of transverse field h(t). a = 0.9(dashed line) 
corresponds to the upper one, a = 0.5(dotted line) corresponds to the middle one, and a = 0.3(solid 
line) is represented by the lower one. We take N = 500, 8 = 0.0001 and r = 250. Inset shows the 
variation of concurrence as a function of transverse field h(t). 
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Fig. 12. Concurrence and discord are plotted as a function of a. Here decoherence factor D(t) = 
0.7025 was taken from transverse Ising quenching data, r = 250, 8 = 0.0001 and t = 251, i.e., after 
crossing the first Ising quantum critical point at h(t) = 1. 

for transverse quenching across Ising critical points achieved by varying h = 1 — t/r. 
Fig. (fT0| shows that discord profile has repeated collapses and revivals similar to 
that of decoherence factor as described in the paper by Damski et al SHI. We see 
complete revival in the region bounded by the Ising critical points h = ±1, while 
beyond the critical point at h = — 1, one observes partial revivals. It is important 
to note that the discord between the two central spins does never vanish even if we 
start with an initial state with a — > 0, when the qubits are in a complete mixed state. 
The qualitative profile of concurrence as a function of h(t) is shown in the inset; 
this turns out to be similar to that of discord. The quantitative difference when 
compared to the behavior of discord is the followoing: concurrence has a higher 
magnitude compared to that of discord for the same value of a and secondly, the 
concurrence between two spins vanishes completely after crossing the second QCP 
after certain value of a which is unlikely as compared to the discord. Concurrence 
stays zero throughout the quenching time below a certain value of a while discord is 
always non-zero. This is identical to the situation when the concurrence and discord 
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are studied in the final state following a quantum quench discussed in previous 
sections. 

Let us now proceed to the small 5 limit when one finds that both discord and 
concurrence exhibit decay with time after crossing QCP h = 1 which is again similar 
to the behavior of the decoherence factor in the same limit SSI ( see pjg_ (pi}). The 
quantitative difference between two measures can be understood using similar argu- 
ments as given above. For example, in Fig. (|12[) . we compare the behavior of discord 
and concurrence as a function of a after crossing the first QCP, i.e., for h(t) < 1. We 
find that the concurrence becomes linear in a beyond a threshold value of a when 
the central spins have a non-zero concurrence; discord is always non-zero for all 
values of a. These results bear close similarity with the corresponding equilibrium 
study 3H. Identical results are obtained when the environment is quenched across 
anisotropic critical points (ACPs) and multicritical points (MCPs) using appropri- 
ate quenching schemes and interaction Hamiltonians between the qubits and the 
spin chain. 

7. Conclusions 

This review attempts to provide a bridge between non-equilibrium dynamics of 
quantum critical systems and quantum information theory. We have shown that 
quenching through QCPs leads to defects in the final states; these defects in fact 
lead to quantum correlations which are manifested in measures like concurrence and 
discord. At the same time both these measures appear to satisfy the same universal 
scaling relation as that of the defect density. Moreover, under certain conditions 
the final state has zero concurrence but non-zero discord implying that correlation 
generated through quench is entirely quantum mechanical. We believe these obser- 
vations will lead to future studies involving dynamics and quantum information. 

Let us conclude with some generic comments: models discussed in the review 
are integrable and also reducible to decoupled two-level systems for which the non- 
adiabatic transition probabilities can be analytically obtained using the Landau- 
Zener transition formula. Moreover, the reduced matrix has a simplified structure 
due to different symmetry properties of the underlying Hamiltonian. Question there- 
fore remains what would happen when one studies quenching of a non-integrable 
model and probe correlations in the final state. Do concurrence or discord present 
in the final state satisfy a universal scaling relation? If so, is the scaling identical to 
that of defect density? How do they evolve when the system thermalizes following 
a quench? These open questions are to be addressed in future studies in this fertile 
area of research. 
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